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Abstract 

We consider a normal-superconducting junction in order to investigate the effect of new physical ingredients on waiting 
times. First, we study the interplay between Andreev and specular scattering at the interface on the distribution of 
waiting times of electrons or holes separately. In that case the distribution is not altered dramatically compared to the 
case of a single quantum channel with a quantum point contact since the interface acts as an Andreev mirror for holes. 
We then consider a fully entangled state originating from spliting of Cooper pairs at the interface and demonstrate a 
significant enhancement of the probability to detect two consecutive electrons in a short time interval. Finally, we discuss 
the electronic waiting time distribution in the more realistic situation of partial entanglement. 
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1. Introduction 

Markus Biittiker was certainly one of the most influen¬ 
tial scientists in the field of mesoscopic physics. Among all 
his important contributions, time in quantum mechanics 
has a peculiar flavor since it occupied his mind at the right 
beginning and at the end of his carrier. Intrigued at first 
by the traversal time of an electron through a tunnel bar¬ 
rier [l|, Q , he came back to this topic after the emer gen ce 
of “on-demand single electron sources” 



IJ , which he greatly contributed to develop 
17, [ll, [I^ [ 2 ^ 1^ . [ 22 . [ 2 ^ via the concept of 
waiting time distribution (WTD) (^, 25. 26|. 

Charge transport at the nanoscale is known to be stochas¬ 
tic due to the quantum nature of particles [16|. There¬ 


fore, going beyond the knowledge of average quantities, 
such as the average electronic current, appears to be un¬ 
avoidable and extremely fruitful at the same time. A deep 
physical insight can indeed be inferred from the fluctua¬ 
tions of the signal and extracted from various observables. 
Noise [l^ and Full Counting Statistics (FCS) [27, 28, 


namely the second moment of current fluctuations and the 
statistics of charges transferred during a long time inter¬ 
val, are among the most popular quantities and have been 


proved to be powerful tools. With the development of 
electron quantum optics and the progress in single 


electron detection at high frequencies 0, 31, 132|, 133|J, it 
is now relevant and possible to consider electron dynam¬ 
ics and time resolved quantities at quantum mechanical 
time scales (typically nano-seconds and below). There¬ 
fore, new theoretical tools have been developed to describe 
the current fluctuations at such time scales, such as fi¬ 
nite frequenc 


37 


cy noise [J, Lka L^, biJ) 1^ 
0, 41, 43 . Wigner functions [4 


36| and FCS 
or the WTD 


The latter. 


describes the statistical distribution of time intervals be¬ 
tween the detection of two electrons and therefore gives ac¬ 
curate information about correlations between subsequent 
electrons. 

The WTD has been studied for particularly simple sys¬ 
tems like single and multiple electronic quantum channels 
connected to two normal leads via a Quantum Point Con¬ 
tact (QPC) [i^, 113 1^. a quantum capacitor [illHi, a 
double quantum dot |45l. 1^ . a train of Lorentzian pulses 
(^ . 1^ or a quantum dot connected to a normal and a 


Email address: inathias.albert@inln.ciirs.fr (M. Albert) 


superconducting lead [4^,1^, among others. In this pa¬ 
per we revisit the physics of Normal-Superconducting (NS) 
junction through the point of view of waiting times in or¬ 
der to illustrate the effect of superconducting correlations 
and entanglement M, M, [Hol. l6Cll| on their distribution. 
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Figure 1: (Color online) Left: Schematic picture of a Normal- 
Superconducting junction. A hole approaching the interface from 
the normal part is either normally reflected or Andreev reflected 
back as an electron. A single electron detector is positioned to de¬ 
tect electrons from Andreev events. Right: energy diagram of the 
setup. The superconducting chemical potential fis is set to an en¬ 
ergy eV above the Fermi energy of the normal part and the gap A 
is much larger than the potential difference eV. 


Indeed, as we will discuss later, such a system may emit 
entangled electrons in the normal part, and leads to inter¬ 
esting features in the WTD. 

The paper is organized as follows. In Sec. [U we de¬ 
scribe the model used for the NS junction and the formal¬ 
ism needed for computing the WTD. In Sec. O we discuss 
the effect of the transparency of the barrier (the energy 
dependence of the Andreev reflection) when the detection 
process is sensitive to only one electronic spin species and 
a certain range of energy. Section |4] is devoted to the effect 
of entanglement between spin up and spin down electrons 
emitted from the superconducting part, on the WTDs. 
We finally conclude and discuss some perspectives in Sec. 
O Moreover, for the sake of clarity, technical details are 
moved to the appendices. Appendix A| demonstrates the 
formal analogy between our setup and a single quantum 
channel conductor for a specific detection process whereas 
important steps for the numerical and analytical calcula¬ 
tions of the WTD in the entangled case are explained in 
[Appendix B[ 


2. Model 

One very important consequence of superconductivity 
is the existence of Andreev reflection. Such a phenomenon 
arises because the superconducting device cannot accom¬ 
modate any single particle excitation with energy below 
the gap A. Therefore, if a single particle like an electron 
or a hole flows from the normal part to the superconduct¬ 
ing part with an energy below this threshold it can only 
be scattered back at the interface. However, there are 
now two possibilities. An electron (a hole) can be either 
normally reflected (specular reflection), that is to say, re¬ 
flected as an electron (a hole) or converted to a hole (an 
electron). This is the so called Andreev reflection which 
originates from the fact that the incoming electron finds 
a partner to create a Cooper pair which can enter in the 
superconductor and leave a hole behind. 


To be more specihc, the system of interest is a polar¬ 
ized NS junction (with an s-wave superconductor), at zero 
temperature, as presented on Fig. [T] The superconduc¬ 
tor chemical potential is set to be at a potential eV 
above the Fermi level Ep of the normal metal. In such 
a situation, there is an incident hole, coming from the 
metal, that can be either normally reflected or Andreev 
reflected as an electron. Another way of picturing the 
Andreev effect is to think about the inverse configuration 
where a Cooper pair in the superconductor (at energy fj,s 
and zero momentum for an s-wave superconductor) splits 
at the interface and gives birth to an entangled pair of 
electrons. From now on, we will take eV much smaller 
than the superconducting gap A in order to focus on this 
sub-gap phenomenon. This also has the benefit to make 
the Andreev time Ia = h/A (the typical time needed for 
an Andreev event) much smaller than t = h/{eV) (the 
typical time separation of two sin gle particle wave pack¬ 
ets emitted in the normal metal 


6l| ). This allows 


us to assume that Andreev events are instantaneous and 
make use of scattering theory. In addition, this assump¬ 
tion allows one to linearize the dispersion relation around 
fis as E{k) = hvpk, with E and k measured from fj,s and 
its corresponding momentum (or around the Fermi level 
since eV -C Ep and ^s)- 

At the interface, the scattering is in general not per¬ 
fect and both normal and Andreev reflection will play a 
role. In order to describe this effect, we use the standard 
Blonder-Tinkham-Klapwijk (BTK) model which has 
been widely used in the literature. The junction is mod¬ 
eled by a point-like barrier potential U(x) = 2ZEpXpS{x), 
where Ai:’ is the Fermi wavelength and Z is a parameter 
measuring the strength of the barrier. It is then possible 
to compute the scattering matrix of this setup exactly and 
obtain the normal and Andreev transmission/reflection co¬ 
efficients 62, 6^. We do not reproduce these results in the 
present paper but give the corresponding numerical values 
of the coefficients when necessary. 

Figure [T] illustrates the scattering processes that we are 
now going to describe mathematically. The incident holes 
of energies fis — E lying between Ep and /i^, arriving from 
the left and propagating to the right will be either normally 
reflected as holes of the same energies with amplitude rjv 
or Andreev reflected as electrons of energies fj,s + E with 
amplitude ta- The incoming scattering state is therefore 
a Slater determinant of holes of the form 57, 58, 591 


eF 

IV'in) = , (I) 

E=0 

where |0) stands for the filled Fermi sea up to fip in the 
normal part. However, in the electron language, this state 
is just the Fermi sea \F) filled up to Ep instead of fig- 
In the following, we will rather use the electronic picture 
to simplify the notation but both pictures are equivalent 
[Hij l . Due to scattering at the interface, the outgoing state 
is therefore a superposition of reflected holes, entangled 
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electrons and non-entangled electrons 56. 57. ^ ^ 62. 


ey 

IV'out) = n (^Ar(^) + 

£;=o 1^) 

X {riy{E) — ■ 

Indeed, it is pretty straightforward to see that the pre¬ 
vious equation, for a given energy, gives birth to three 
kinds of term with different levels of complexity. The terms 
|F) correspond to 

non-entangled contributions whereas ~ 

''1(B) 4.''-fc(-£;) describes fully entangled electrons 

originating from the splitting of a Cooper pair at the in¬ 
terface. When Andreev reflection is absent {va = 0), the 
Fermi Sea is unperturbed by the interface and nothing in¬ 
teresting happens. Counter-intuitively, perfect Andreev 
reflection does not lead to perfect entanglement. On the 
contrary, the state is a Slater determinant of non-entangled 
electrons and the NS junction acts as a conventionnal elec¬ 
tron source [H^. It appears that the maximally entan¬ 
gled situation arises when Andreev and normal reflection 
probabilities are both equal to one half. Nevertheless, the 
WTD of a fully entangled state has never been studied to 
our knowledge and we will take the opportunity to study 
it in this paper before considering the general and more 
realistic state emitted at the interface. 

In order to conclude this section, we recall a few defini¬ 
tions about WTDs. As mentioned in the introduction, the 
waiting time r is defined as the time delay between the 
detection of two single particles. Due to scattering and 
the quantum nature of particles, this time is a random 
variable, which distribution (the WTD) brings an elegant 
and instructive picture of the physics. For stationary sys¬ 
tems, namely when there is no explicit time dependence, 
the WTD yV(r) depends on r only (and not on absolute 
time) and is closely related to the Idle Time Probability 
(ITP) n(T), the probability to detect no electron during a 
time interval r 


W(r) = (r) 


d^n(r) 

dr 2 


( 3 ) 


dn(r) 

dr 


-1 


T=0 


is the mean waiting time |4S 


where (r) = — 

To go further, we must now specify the detection proce¬ 
dure to compute the WTD. In what follows we will as¬ 
sume perfect single electron projective measurement but 
will consider two different situations as described below. 
Such perfect detection process is still theoretical and ex¬ 
tremely challenging at very short time scales (nano-second 
and below) but recent experiments [i, 32, ^ are very 

promising about this issue. In order to compute the ITP, 
we regularize the scattering problem as usual . We dis¬ 
cretize energy, ranging from 0 to eV, in N slices and wave 
vectors as kn = ^md consider the limit iV —>• oo to 

mimic a stationary process. 



Figure 2: (Color online) WTD of Andreev events for various barrier 
transparencies. As Z is increased, Ra decreases from one to zero as 
shown in the inset. For Z = 0, Ra = 1 and the WTD is described by 
the Wigner distribution (blue dots, Eq. 0). Close to perfect spec¬ 
ular reflection {Ra "C 1) the distribution approaches an exponential 
law except for very small times (black squares). For intermediate Z, 
small oscillations with period r are superimposed to an exponential 
decay. Together, with the dip at r = 0, they are manifestations of 
Pauli’s exclusion principle 1^ . 


3. Andreev mirror for detecting WTD of holes 


The first application of this setup will be to measure 
the WTD of holes via the detection of Andreev reflected 
electrons or in other words the WTD of Andreev events 
641 . A single electron detector is located at a position xq 
far away from the interface. In this section we assume that 
the detector measures only one spin orientation that we 
will choose upward for concreteness. Moreover, we make 
the additional assumption that it is only sensitive to ener¬ 
gies above using a quantum dot for instance [H^. As 
a consequence, this allows us to avoid complications due 
to entanglement Ig^^which will be the subject of the 
next section. Following , the ITP is defined as 


n^(T) = (^outl : e '5t,E>(^s : \tjjont) (4) 

with the explicit condition that only electrons with ener¬ 
gies above /is contribute. Here Qt = Ct(^) dx 

and : • • • : stands for the normal ordering. In principle, the 
ITP cannot be reduced to a single determinant since the 
many body state is not a Slater determinant. However, 
according to the assumptions introduced above, only one 
term survives and the hnal result can be cast as a single de¬ 
terminant [ 2 ^, 0,1^, [ 4 ^ . Under these conditions, entan¬ 
glement no longer plays a role and the problem thus boils 
down to a single quantum channel with energy dependent 
transmission |49l| . where the role of the energy-dependent 
transparency t{E) of the QPC is played here by rA{E), the 
Andreev reflection amplitude (see |Appendix A| for details). 

In Fig. [21 we have plotted the spin up electronic WTD 
for different barriers strengths Z = 0 (perfect Andreev 
reflection), Z = 0.5 and Z = 1.5 (strong barrier) as a 
function of t jr. For information we show as an inset the 
corresponding Andreev reflection coefficient Ka = 
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as a function of energy. However, the energy dependence 
is very weak since the energy window eV is supposed to 
be much smaller than the superconducting gap A. As 
expected, the WTD for Z = 0 is approximately given by 
the Wigner surmise [1^ 


32 

Wwd('r) = ^ =Texp 
TT^ T 


.1 (LV 

TT \T/ 


(5) 


Indeed, in that case the train of free holes is perfectly 
converted at the interface into free electrons which are 
described by random matrix theory 25, As Z is 


increased, the distribution is broadened since Ra is no 
longer equal to one and therefore not all holes are con¬ 
verted into electrons. The situation is exactly equivalent 
to free electrons injected into a single quantum channel 
and partitioned by a quantum point contact with trans¬ 
mission probability Ra due to the electron-hole symmetry 
in the system. The difference is only conceptual since here 
the detector is measuring indirectly the statistics of holes 
converted into electrons by the NS junction that acts as an 
Andreev mirror. Finally, for large Z, namely, small An¬ 
dreev reflection, most of the holes are normally reflected 
and the detector collects rare events which are almost un¬ 
correlated and the WTD is exponential (except for very 
short times). Indeed, following we derive the asymp¬ 
totic behavior of the WTD in the long time limit. For 
Ra = 1 the decay is Gaussian with algebraic corrections 
well described by the Wigner surmise and for partial An¬ 
dreev reflection {Ra < 1) it is exponential with a rate 
given by the geometrical mean of the logarithm of (1 — Ra) 
over energy in [^s, ^J-S + eH], namely 


W(r) 


exp 


ln(l - Ra) r/r 


gi^/r), 


( 6 ) 


where g{y) is an oscillatory function that decays as l/y"^ 
and depends on In^ (1 — Ra ). 

To conclude this section, we note that the role of elec¬ 
trons and holes may be interchanged by inverting the po¬ 
larization (eV —> —eV) due to electron-hole symmetry. 


4. Effect of entanglement on waiting times 

We now move a step forward and discuss the effect 
of entanglement on waiting times. As mentioned before, 
Andreev reflection might be thought in terms of splitting 
of a Cooper pairs in the vicinity of the interface, leading 
to the injection of two entangled electrons with opposite 
spins and energy (with respect to ys) in the normal part 
as shown in Fig. [3] (left part). However, we have seen that 
the scattering state (Eq. ([3])) is a mixture of entangled 
and non-entangled components which makes the effect of 
entanglement hard to separate from other physical ingre¬ 
dients. Therefore, we focus here on the WTD of the fully 
entangled state of the form 


Entangled 



Non — Entangled 



Figure 3: Schematic picture of the setup. Left: a Cooper pair is split 
at the interface of a NS junction and gives birth to two entangled 
electrons in the normal part. Right: two independent electrons as 
they would be emitted if the right part was also a normal metal. 


eV 


IV'fe) = R 


E=0 


J J 


- d d 


V2 


\F)^ 


^ (7) 

even if this is not the real quantum state emitted at the 
interface. Indeed, this work on NS junction rises a funda¬ 
mental question of the effect of entanglement on waiting 
times that has never been considered to our knowledge 
and deserves to be investigated with proper care. At the 
end of this section we will give a few hints on how the 
additional terms included in the full state (Eq. ([2])) mod¬ 
ify the picture but the common thread of this section will 
be the study of the fully entangled quantum state (Eq. 
©)■ To point out the specific features of entanglement, 
we will compare this situation to the one of two indepen¬ 
dent electrons which would correspond to the emission of 
two electrons with opposite spins from a normal metal (see 
Fig. [3] right) [i^. This time, the single electron detector 
is sensitive to both spins and energies between /is — eV 
and /is -|- eV. Using the energy discretization mentioned 
before, the entangled state reads 


N 


IV'fe) = R 


n=l 


j j 

-fend 


— c 




4.^-fcn.t 


V2 


\F ), (8) 


whereas non-entangled electrons emitted from a normal 
lead in an energy window eV (see Fig. [3] right) would be 
described by IV'ne) = OLi 

According to the detection process mentioned above, 
the ITP is formally given by [2g 


n(T) = (: ; e-'Si y (g) 

where the quantum average is taken over the state Id^FE) or 
I^Tne) leading to Hfe and Hne respectively. Without any 
further assumption it is clear that the final result cannot 
be expressed as a single determinant but rather as a sum 
of 2^^ terms except in the non-entangled case where the 
ITP factorizes to 
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HneM = (: e-«t :) (: 


3 Q-i- A - 


:) = det(l - Qt) det(l - Qi), 

( 10 ) 

where the averages are taken over the two spin sectors 
separately [i^. Indeed, each term in the product of Eq. 


can be split into two parts: A\^ = 


V2 


and 


Bi = -- 


-■ IV'fe) is then a sum of 2^ terms of the 


V2 

form Oz^iCilE), where C„ can be either An or Bn- We 
shall denote a particular term by a string made of a suc¬ 
cession of A’s and B’s defined as follows: |A, B, = 

j^^B\...B^j^\F). This is a Slater determinant, while I^Te) 
can not generally be cast as a simple Slater determinant. 
IIfe will then be the sum of 2^^ terms of the form, typ¬ 
ically, Tij = ;; with 

and two generally different Slater determinants. As 

previously shown in Refs. , each Tij is also a de¬ 

terminant. Details about the procedure to calculate the 
ITP are given in [Appendix B[ 


4-1- Numerical results for the fully entangled case 

We now discuss our results obtained from a direct and 
exact enumeration of the 2^^ terms in the ITP. Owing 
to the exponentially growing number of terms, this ap¬ 
proach is limited to relatively small values of N. Figure 
0^) presents the WTD for increasing values of N up to 
N = 12. As can be seen, the curves reasonably con¬ 
verge to a limiting distribution that would be obtained 
for N —> -|-oo. To insure this, we have computed several 
finite size corrections that we will discuss later on. Figure 
|3 )d) compares the WTDs of entangled and non-entangled 
electrons which are qualitatively similar. However, as 
we will discuss in more detail in the next subsection, the 
maximum of the curve in the entangled case is more pro¬ 
nounced and closer to zero than in the independent case. 
The presence of such a peak in the WTD is the hallmark 
of pair rigidity due to entanglement. To be more quantita¬ 
tive we can evaluate the probability that the waiting time 
is smaller than the average waiting time h/ 2 eF (which is 
the same in both situations). We hnd that this probability 
is about ten percent larger in the entangled case, demon¬ 
strating that the entangled electrons are more correlated 
than the non-entangled ones. 


Jf..2. Short time behavior for the fully entangled case 

The short time behavior is one of the most interesting 
characteristics of the WTD as it reflects the short time cor¬ 
relations encoded in the many body state and not the ones 
due to scattering. For free electrons it is universal since it 
is the expression of the Pauli’s principle. Indeed, two spin¬ 
less electrons cannot be emitted in the same state which 
enforces the WTD to start from zero (with a quadratic 
behavior). However, if the electrons have other degrees of 
freedom like spin or if the mesoscopic conductor supports 
several channels, the WTD may start from a non-zero 



Figure 4: (Color online) WTD of fully entangled electrons, a) WTDs 
for increasing values of N. The inset shows how VV(0) converges to 
the asymptotic value 1 /t as predicted by Eq. 1111 (full black line), 
b) Comparison between the WTD of entangled electrons and non- 
entangled electrons (see text). Inset: short time behavior in the 
entangled case for N = 12 (full red line and green dashed line (see 
Eq. Him and N —)■ +oo (black long dashed lines). 


value but correlations are still visible and universal (d^ . 
However, these correlations only originate from electrons 
of the same channel and not between different channels. 
They eventually disap pea r in the limit of large number of 
independent channels |49l| . 

In order to get the short time expansion of the WTD 
up to second order in time, we need to expand the ITP 
(Eq. (O) to fourth order in terms of moments of Q 4 ,, 
and their products. This is a straightforward but some¬ 
how cumbersome calculation that we do not reproduce in 
detail here. After some algebra we obtain the short time 
expansion of the WTD, including finite size corrections. 
The hnal expression for t <C t reads, to second order in r 


W(t) 


1 

r 






( 11 ) 


where Cn = 4(1-1-^) in the entangled case and one in the 
absence of NS junction. The insets of FigH^) and Figl 4 l 3 ) 
show how this prediction is conhrmed by our numerical 
evaluation of the WTD in the entangled case. In both 
situations, the WTD starts from the same initial value 
(in the N ^ 00 limit) with a quadratic behavior but the 
coefficient is four times larger in the entangled case. Again, 
this reflects the pair rigidity due to entanglement. 
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4-3. Long time limit for the fully entangled case 


Differences between entangled and non-entangled elec¬ 
trons can also be inferred from the tail of the WTD. Fol¬ 
lowing the approach developed in ^ we obtain the asymp¬ 
totic behaviors of the WTD in the limit r ^ r (for details 
see [Appendix B[). In both cases we find a Gaussian decay 


n(r) 


exp 



( 12 ) 


with c a constant which is equal to one in the entangled 
case and two in the non-entangled case. Therefore, the 
ITP or equivalently the WTD decays more slowly in the 
entangled case, meaning that the two electrons are not at 
all independent. 


4 .4- Realistic situation in the tunneling limit 

As mentionned before, entangled pairs are only one 
component of the real out-going scattering state which ap¬ 
pears to be much more complex. For abritrary value of the 
reflection coefficients Rn and Ra, numerical evaluation of 
the WTD is pretty challenging since the number of terms 
grows as In that case, direct numerical evaluations 
are restricted to very small values of N (typically 4 or 5) 
which is not sufficient to mimic a stationary situation. It is 
then necessary to resort to more sophisticated approaches 
like Monte-Carlo sampling of the ITP for instance. How¬ 
ever, there is a simple but rather interesting limit which 
is amenable to analytical calculations, namely the tunnel¬ 
ing limit for Ra <C 1 (i?Ar ~ G. In that case, the full 
many-body state simplifies to |58l. l56| 


l^t) = 


eV 


1 + rain 


E=0 


\F), (13) 


namely, to first non-trivial order in Ra , it is made of scarce 
singlet pairs. Following the same approach as before we 
obtain the small time behavior 


5. Conclusion 

We have presented a theory of waiting times in polar¬ 
ized NS junction at zero temperature. In this setup, An¬ 
dreev reflection brings new characteristic features in the 
WTD such as entanglement. 

If a detector is sensitive to only one type of spin and 
to energies above the superconductor chemical potential, 
the situation is reminiscent of a single quantum channel 
connected to two normal leads via a QPC with energy- 
dependent transmission. The interface acts as an Andreev 
mirror and allows to measure the WTD of holes converted 
to electrons. However, if the detector measures electrons 
with both spins in the whole energy window above the 
Fermi sea of the normal conductor, entanglement between 
electrons leaves fingerprints in the WTD. Although still 
academical, we have focused on a fully entangled state 
which is only one component of the full many-body state 
flowing out of the NS interface. In that case, such signa¬ 
tures are visible for both small and large waiting times but 
the most important feature is the existence of a peak in 
the WTD centered before the average waiting time. When 
taking into account all components of the scattering state, 
we have shown that some charateristics are still visible in 
the WTD in the tunneling limit. 

In the near future, it would then be useful to extend 
this work beyond the ideal situation and the tunneling 
limit and evaluate the effect of entanglement for arbitrary 
values of reflection coefficients. In addition, a natural ex¬ 
tensions would be to study correlations between waiting 
times of different spin species in the spirit of (H^ . This 
would probably yield an even clearer signature of entan¬ 
glement than the WTD itself. 

Among other future investigations, it would be pos¬ 
sible to study the effect of cross-Andreev reflection in a 
Superconducting-Normal-Superconducting junction or the 
physics brought by exotic states like Majorana modes cre¬ 
ated by Majorana guns 


W(t) 




and the long time asymptotics 


W(r) ^ exp[—i?AT/'r]. 


(14) 


(15) 


At short waiting times, the distribution starts from a 
constant value and the physics is dominated by correla¬ 
tions within a single Cooper pair. Indeed, simple consid¬ 
erations show that pairs are roughly separated in time by 
T/Ra whereas electrons from the same pair are rather in 
a span of time of the order of r. The situation is then 
akin to a single pair which has been studied by Hassler 
et al. [ 3 ^ and confirms our predictions. At large time, 
the decay is exponential with a rate twice smaller than in 
the non-entangled case, which again is an hallmark of pair 
rigidity. 
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Appendix A. Supplementay information on the An¬ 
dreev mirror 

This appendix is dedicated to prove the formal corre¬ 
spondence between the WTD of a single channel normal 
conductor and the Andreev mirror defined in section [3] 
This happens because the detection process is only sensi¬ 
tive to spin up electrons and energies above /ig, the su¬ 
perconducting chemical potential. To begin the proof we 
recall the definition of the ITP in terms of the scattering 
state 


N N 

(-P'1 -f rAcl,^^)\F). 

2=1 7=1 

(A.6) 

We then recognize the expression for the ITP of a spin¬ 
less quantum channel with bias et^ and transmission am¬ 
plitude rA [13 ■ 

Appendix B. Fully entangled case 


nt(r) = (V’outI : (A.l) 

with 

N 

IV'out) = - r*Acl^^^cl^^^^^)\F) 

2 = 1 


(A.2) 


Appendix B.l. Numerical procedure in the fully entangled 
case 

In this Appendix we give some details on the construc¬ 
tion on n(r). Starting from the N-body states reads 


N 


w = n 


n—l 


(cL,tcLfe„q 




V 2 


\F), (B.l) 


where N corresponds to the number of slices due to the 
energy discretization between p,s and fis + eF (which has 
to be taken equal to infinity at the end of the calculation) 
and k is measured with respect to ks the Fermi momen¬ 
tum at energy fis- There are two important steps in the 
derivation. The first one consists in expanding the prod¬ 
ucts and rearranging the terms in order make the relevant 
terms appear (the ones with spin up and positive energy 
electrons). If we restrict ourselves to these terms we build 
terms like \C) = ni!=i(^v + ■ Then, we 

include the other terms of the product and show that they 
are irrelevant. Indeed, in the following example 


(C|(r)V - rAC-k^,tCkui) ■■ e : (r^ - j^cLfe^_|)|C), 

(A.3) 

the two off-diagonal terms vanish since the operator : 
conserves the number of particles. Concerning the two di¬ 
agonal ones, they leave us with 


(^ATp + : e : \C) = (C| : e : \C). 

(A.4) 

Then, we can repeat the same argumentation for all the 
other irrelevant terms from z = 1 to A. Finally, we are 
left with 


N N 

{F\\{v\-.e-^^A>^s ■.^V,\F) (A.5) 

i=l 


with \F) is the Fermi sea defined in the main text and 


27rn |eF| 
N hvp 


(B.2) 


The probability to detect nothing in a range of time r with 
a detector sensitive to all the energies and all spins can be 
written as 


n(T) = (: e-^^ :: :), (B.3) 


where the average is taken over lip). From here, it is con¬ 


venient to separate Eq. (EH) in two parts cj. 






.,t -fend 


and 


t'^-fen 




in order to get a sum of 2^ terms where each 
term is a Slater determinant. To simplify, it is useful to 
change the notation; now each configuration is associated 
to a specific ket (bra) and each element of this configu¬ 
ration can be mapped to an Ising classical spin. Namely, 
ct I (ct j .) is associated to ai = -|-1(—1). Eq. 




• (IB 3p becomes [25l. l3 


n(^) = ^ E {WUM\{a}k). (B.4) 

W}b,W}k 


where the summation runs over the 2 ^^ configurations 
for the bra/ket (One example of possible configuration is 
((TifT 2 (T 3 | = (1, —1,1| for a system consisting of three elec- 
tron pairs). Each element of the matrix M can be written 
as [39| 

_ r X/Af) j. V/AT\ fT3 

^n,m 


with Vi = {r’^ + rAc],, ^)- As a second step, we with = cr„n — and X = r/f. 
rearrange the Eock space by moving all the spin f to the 
left. Spin I are invariant under : e“‘5t,fc>fcs :, thus we only 
have to compute 
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Appendix B.2. Calculation o/n(r) in the long time limit 
We now explain how to compute the long time asymp¬ 
totics of the ITP of the fully entangled state. In that case, 
the matrix reduces to usual Toeplitz matrix which inter¬ 
venes in the statistics of levels in the Gaussian Unitary 
Ensemble ensemble of random matrices [69| . Namely, 

M = I — exp{iipn) sm{(pn)/nTr, (B.6) 

with tpn = mrX/N and X = t jr (except for an unim¬ 
portant factor TT in the angle Pn)- The asymptotic behav¬ 
ior of the WTD for large r is equivalent to the statistics 
of having no level in a large energy range. Using this 
correspondence, the main terms read, for large W(r) ~ 
exp[—(r/r)^/2], with C a rather complicated 
constant which can be find in reference [z3- The decay is 
Gaussian with algebraic corrections. 
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